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Abstract 



' Addition theorems can be constructed by doing three-dimensional Taylor expansions according 

, to /(r + r') = exp(r' • V)/(r). Since, however, one is normally interested in addition theorems 

Oh' of irreducible spherical tensors, the application of the translation operator in its Cartesian form 

fH , exp{x'd/dx)exp{y'd/dy)exp{z'd/dz) would lead to enormous technical problems. A better alter- 

native consists in using a series expansion for the translation operator exp(r' • V) involving powers 
of the Laplacian and spherical tensor gradient operators 3^™(V), which are irreducible spherical 
tensors of ranks zero and £, respectively [F.D. Santos, Nucl. Phys. A 212, 341 (1973)]. In this way, 
it is indeed possible to derive addition theorems by doing three-dimensional Taylor expansions [E.J. 
Weniger, Int. J. Quantum Chem. 76, 280 (2000)]. The application of the translation operator in 
its spherical form is particularly simple in the case of B functions and leads to an addition theorem 
' with a comparatively compact structure. Since other exponentially decaying functions like Slater- 

type functions, bound-state hydrogenic eigenfunctions, and other functions based on generalized 
Laguerre polynomials can be expressed by simple finite sums of B functions, the addition theorems 
for these functions can be written down immediately. 

1 Introduction with r, r' G in terms of other functions that 

only depend on either r or r'. 
In many branches of physics and physical chem- atomic and molecular calculations, one is 

istry - for example in electrodynamics Q, in das- usually interested in irreducible spherical tensors 
sical field theory @ , or in the theory of intermolec- ^^lich can be partitioned into a radial part and an 
ular forces § - an essential step towards a solution angular part that is given by a spherical harmonic: 



Ftir) = h{r)Yr{r/r). (1.1) 



of the problem under consideration consists in a 
separation of the variables. 

Principal tools, which can accomplish such Thus, the function /(r it r'), which is to be ex- 
a separation of variables, are addition theorems: panded, as well as the functions, that only depend 
These are expansions of a given function /(ritr') on either r or r', should be irreducible spherical 
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tensors. 

The best known example of such an addition 
theorem is the Laplace expansion of the Coulomb 
potential in terms of spherical harmonics: 



1 



E 

l=Q m=-t 



47r 



2^ + 1 



[yr(r</r<)]*yr(r>/r>), 



r< = min(r, r') , r> = max(r, r') . 



il.2) 



The Laplace expansion leads to a separation of 
the variables r and r'. However, the right-hand 
side of this expansion depends on r and r' only 
indirectly via the vectors r< and r> which satisfy 
|r<[ < |r>[. Hence, the Laplace expansion has 
a two-range form, depending on the relative size 
of r and r'. This is a complication which occurs 
frequently among addition theorems. 

There is an extensive literature on addition 
theorems. Since addition theorems can be viewed 
as expansions in terms of spherical harmonics with 
arguments r/r and r'/r', they are often treated in 
books on angular momentum theory or on related 
topics (see for example pp. 163 - 169 of Q or Ap- 
pendix H.4 of II). 

Particularly well studied are the addition theo- 
rems of the solutions of the homogeneous Laplace 



equation 



6^ 



21[1, of the homogeneous Helmholtz 
and of the homogeneous 



equation [|r3,|l|,||-|| 
modified Helmholtz equation [lOt 26|-28,|3C|. This 
large number of references does not only reflect the 
importance of these functions, but also the rela- 
tive ease with which their addition theorems can 
be derived. The derivation of addition theorems 
for other functions, which are not solutions of the 
equations listed above, is according to experience 
significantly more difficult. 

In the context of atomic and molecular elec- 
tronic structure calculations, addition theorems 
have traditionally been used for the separation 
of the variables of interelectronic repulsion inte- 
grals and of other multicenter integrals. For that 
purpose, one also needs addition theorems for the 
so-called basis functions. Probably the oldest at- 
tempt is the zeta function method of Barnett and 



Coulson 1 31 1 which tries to obtain addition theo- 
rems for exponentially decaying functions by ap- 
plying suitable generating differential operators to 
the well-known addition theorem of the Yukawa 
potential e~"^/r. In the case of Slater-type func- 
tions, this approach did not lead to a complete 
success: It was not possible to do the differentia- 
tions in closed form, and the coefficients occurring 
in the zeta function expansion had to be computed 



recursively 1 31 , 32 ] . 

Another general method for the derivation of 
addition theorems is the alpha function method, 
which was introduced by Per-Olof Lowdin in his 
seminal paper on the quantum theory of cohesive 
properties of solids [^3|, and which was later used 
and extended by numerous other authors |34-]50[. 

In 1967, Ruedenberg |51] and Silverstone [^2| 
suggested independently to derive addition theo- 
rems via Fourier transformation. If the function 
/(r lb r') is an irreducible spherical tensor of the 
type of (|l.lD , then the angular integrations can 
be done in closed form, and there only remain 
some radial integrals in momentum space contain- 
ing two spherical Bessel functions (see for exam- 
ple Eq. (7.5) of (5^). Unfortunately, these inte- 
grals can be very difficult. There is the additional 
complication that even if explicit expressions for 
these integrals can be found, the radial variables 
r and r' need not be separated in the final result. 
Thus, the success of this approach depends cru- 
cially upon one's ability of handling complicated 
integrals containing products of spherical Bessel 
functions. 

Addition theorems can be defined as three- 
dimensional Taylor expansions (see for example p. 
181 of i5l): 



/(r + rO = f;^^^/(r) 



= e 



n=0 
r'-V 



n! 
/(r) 



(1.3) 



Thus, the translation operator e"" generates 
f{r + r') by doing a three-dimensional Taylor ex- 
pansion of / around r. 

From a practical point of view, three- 
dimensional Taylor expansions do not seem to be 
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very useful for the derivation of addition theorems. 
In atomic and molecular physics, one is usually in- 
terested in addition theorems of irreducible spher- 
ical tensors of the type of (1.1) which are defined 
in terms of spherical polar coordinates r, 9, and 
(f). Accordingly, differentiations with respect to 
the Cartesian components x, y, and z of r would 
lead to extremely messy expressions. Moreover, 
it is highly desirable to express the angular part 
of the addition theorem in terms of the spheri- 
cal harmonics Yj^^(r/r) and Y^^^(r'/r'). Thus, 
any attempt to derive an addition theorem for 
an irreducible spherical tensor via a straightfor- 
ward application of the translation operator in its 
Cartesian form e^'d/dx^y'd/dy^z'd/dz ^^^j^ ^g^^j 

enormous technical problems. 

Nevertheless, it is not only possible to obtain 
addition theorems via three-dimensional Taylor 
expansions, but this approach may even be the 
most natural and the most widely applicable tech- 
nique . This is accomplished by expanding the 
translation operator e"" in terms of differential 
operators that are irreducible spherical tensors. 
Then, the angular part of the addition theorem 
is obtained via angular momentum coupling, and 
only radial differentiations have to be done explic- 
itly. Obviously, this leads to a significant simplifi- 
cation p5|]. 

In [p5|, it was shown how the Laplace expan- 
sion (|l]^) of the Coulomb potential and the ad- 
dition theorems of the regular and irregular solid 
harmonics and of the Yukawa potential e~°'^ /r can 
be derived by applying the translation operator in 
its spherical form. Of course, these addition the- 
orems are all comparatively simple, and they do 
not necessarily provide convincing evidence that 
the Taylor expansion method is indeed a useful 
tool for the derivation of more complicated addi- 
tion theorems. However, on p. 293 of ||5^, is was 
emphasized that this method had been applied 
successfully also in the case of other functions of 
physical interest. Thus, it is the intention of this 
article to provide some evidence that these claims 
are indeed true (further articles on additions the- 
orems are in preparation |5^). Moreover, in 
the addition theorem of the so-called spherical La- 



guerre Gaussian functions was derived by applying 



the techniques described in |55|. 

Currently, the vast majority of molecular elec- 
tronic structure calculations use Gaussian basis 
functions although these functions are neither able 
to reproduce the cusps |58| of the exact solu- 



tions of molecular Schrodinger equations at the 
nuclei nor their exponential decay [5^, 6C | . The 
only, but nevertheless decisive advantage of Gaus- 
sian functions has been the relative ease, with 
which the molecular multicenter integrals can be 
computed. Nevertheless, there are still many re- 
searchers who hope that ultimately the unphys- 
ical Gaussian functions can be replaced by the 
physically better motivated exponentially decay- 
ing functions and who are doing research in this 
direction. Accordingly, there is a large number of 
recent articles on molecular multicenter integrals 
of exponentially decaying functions |6l|- ^] . 

So far. Slater-type functions have been the 
most important and most widely used exponen- 
tially decaying basis functions in atomic, molec- 
ular, and solid state calculations. Unfortunately, 
the molecular multicenter integrals of Slater-type 
functions are notoriously difficult, and so is its ad- 
dition theorem which can be used for the separa- 
tion of variables in multicenter integrals. Over the 
years, a large variety of different techniques were 
used for the construction of addition theorems for 
Slater-type functions fS^^, ||, ||, H, H, ||, ||, 
^,9C, 91|, and most likely there is still room for 
improvements. 

The topic of this article is not the derivation 
of addition theorems for Slater-type functions but 
for another class of exponentially decaying func- 
tions, the so-called B functions At first sight, 
this may look surprising since B functions have 
a comparatively complicated mathematical struc- 
ture, and it is by no means obvious that any- 
thing can be gained by considering B functions 
instead of the apparently much simpler Slater- 
type functions. However, B functions have some 
remarkable mathematical properties which give 
them a unique position among exponentially de- 
caying functions and make them especially useful 
for molecular calculations. Due to these ad van- 
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tageous properties it is not only relatively easy 
to derive an addition theorem with the help of 
three-dimensional Taylor expansions as described 
in ||5^, but the addition theorems for B functions 
also have comparatively simple structures [p^. 

Moreover, all the commonly used exponen- 
tially declining functions as for example Slater- 
type functions or bound-state hydrogenic eigen- 
functions can be expressed as simple finite sums 
of B functions (see for example Section III of |^3| 
and references therein). These advantageous fea- 
tures have led to a considerable amount of recent 
research on B functions and their multicenter in- 



tegrals [62,68-31,87 



Addition theorems for B functions have al- 
ready been derived in |Q by applying techniques 
that closely resemble the zeta function method of 
Barnett and Coulson |Q. So, the emphasis of this 
article is not so much on the derivation of new re- 
sults but on the discussion of methodological ques- 
tions related to the application of the translation 
operator in its spherical form to exponentially de- 
caying functions. 

In Section ^, the definitions and conventions 
used in this article are listed. In Section ^ the 
expansion of the translation operator in terms 
of irreducible spherical tensors is reviewed, and 
the properties of the differential operator 3^^(V), 
which is the key quantity in this expansion, are 
discussed. Section ^ covers B functions and their 
mathematical properties, and in Section ^, an ad- 
dition theorem for B functions is derived by doing 
a three-dimensional Taylor expansion. Finally, in 
Section ^ there is a short discussion of the results 
presented in this article. 



2 Definitions 

For the commonly occurring special functions of 
mathematical physics we use the notation of Mag- 
nus, Oberhettinger, and Soni |94| unless explicitly 
stated. 

For the set of positive and negative integers, 
we write Z = {0, ±1, ±2, . . . }, for the set of pos- 
itive integers, we write N = {1, 2, 3, . . . }, and for 



the set of non-negative integers, we write No = 
{0, 1,2,...}. The real numbers are denoted by M, 
and for the set of three-dimensional vectors with 
real components we write M^. 

For the spherical harmonics Yf^{9,(l)), we use 
the phase convention of Condon and Shortley , 
i.e., they are defined by (p. 69 of [E3) 



m+\m\ 



{2i+l){i- \m\)l 
47r(^+ |m|)! 



1/2 



(2.1) 



Here, Pj"^'(cos6') is an associated Legendre poly- 
nomial (p. 155 of |l96|): 

= (l-x2)-/2__P,(^). (2.2) 

Under complex conjugation, the spherical har- 
monics satisfy (p. 69 of jS^): 



(2.3) 



For the regular solid harmonics, which is a so- 
lution of the homogeneous Laplace equation, we 
write: 



(2.4) 



The regular solid harmonic is a homogeneous poly- 
nomial of degree I in the Cartesian components x, 



y, and z of r (p. 71 of [p4[): 
'2£ + l 



yrir) 



An 



{e + m)i{e-my. 



1/2 



E{-x — iy)'^^^{x - iyYz^ m 2fc 



A:>0 



2»n+2fc(^ k)\k\{i - m - 2A;)! 



This formula holds for the Cartesian components 
of an arbitrary three-dimensional vector. Conse- 
quently, it can be used to define the differential 
operator 3^™(V) by replacing in (|2.5| ) the Carte- 
sian components x, y, and z of r by the Cartesian 
components d/dx^ d/dy, and d/dz of V. 
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For the integral of the product of three spher- 
ical harmonics over the surface of the unit sphere 
in M^, the so-called Gaunt coefficient we write 

(£37713 |£2m2|^imi) 

(i^)]*y,7^(i7)y,7(i7)do. (2.6) 

The Gaunt coefficient can be expressed in terms 
of 3jm symbols (see for example p. 168 of 



{i3m3\£2m2\£imi) 



(-1) 



il h 









47r 

mi m2 



-1713 



(2.7) 



It follows from the orthonormality of the spheri- 
cal harmonics that the Gaunt coefficients linearize 
the product of two spherical harmonics: 



e=e, 



^"^^ {imi + m2\£imi\£2m2) 

(0) . (2. 



The symbol X^*-^^ indicates that the summation 
proceeds in steps of two. The summation limits 
in (p.^), which follow from the selection rules sat- 
isfied by the 3jm symbols, are given by p8| 



in 



if ^max + Amin is even, 



if Cax + Amin IS odd, 



(2. 



where 
A, 



max{\li- e2\,mi+m2) . (2.10) 



It may be of interest to note that several arti- 
cles on Gaunt coefficients have appeared recently 



In the sequel, the following abbreviations will 
be used: 

A£ = (£i + £2-^)/2, (2.11) 
Ail = ii-ii+i2)/2, (2.12) 
Ai2 = i£ + ii-i2)/2, (2.13) 
a{i) = {ii+i2 + e)/2. (2.14) 

If the three orbital angular momentum quantum 
numbers ii, £2, and i satisfy the summation lim- 
its ( |2.9D , then these quantities are always positive 
integers or zero. 



3 The Translation Operator in 
Spherical Form 

In this Section, the expansion of the translation 
operator e"" in terms of irreducible spherical ten- 
sors will be reviewed. Moreover, the for our pur- 
poses most important properties of the differen- 
tial operator 3^^(V), which is the key quantity in 
the expansion mentioned above, will be discussed 
shortly. Additional details be found in Sections 3 



and 4 of |55||. 

If we want to derive an addition theorem for 
some function / by applying the translation oper- 
ator according to (|l.3| ) , / has to be analytic at the 
expansion point r as well as at the shifted argu- 
ment r + r'. Consequently, (|l]^) is in general not 
symmetric with respect to an interchange of r and 
r'. Thus, if the three-dimensional Taylor expan- 



sion (|l.3| ) converges, it is in general not possible 
to expand / around r' and use r as the shift vector 
according to 



/(r + r') 



(r-V 



fir') 



n=0 
„r-V' 



fir') 



(3.1) 



PHio^] 



The reason is that (|l.3|) and (|3.lD can only hold 
simultaneously for essentially arbitrary vectors 
r, r' G if / is analytic at r, r', and at r + r'. 
However, most functions, that are of interest in 
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the context of atomic and molecular quantum me- 
chanics, are either singular at the origin or are not 
analytic at the origin. Obvious examples are the 
Coulomb potential, which is singular at the origin, 
or the Is hydrogen eigenfunction, which possesses 
a cusp at the origin In fact, all the commonly 
used exponentially decaying functions as for exam- 
ple Slater-type functions or also B functions are 
not analytic at the origin. 

The reason for the non-analyticity is that the 

1 /2 

three-dimensional distance r = [x^ + y'^ -\- z^] is 
not analytic with respect to x, y, and z at the ori- 
gin r = 0. In contrast, = + + z'^ and the 
regular solid harmonic 3^™(r) are analytic since 
they are polynomial in x, y, and z. Consequently, 
a Is Gaussian function exp(— ar'^) is analytic at 
r = 0, but a Is Slater-type function exp(— ar) is 
not. 

Thus, for the derivation of addition theorems 
for B functions, which are not analytic at the ori- 
gin, we have to use the translation operator in the 
following form. 



/(r<+r>) 



(r< • V>) 
n! 



n=0 

e'-<-^>/(r>). 



(3.2) 



where < |r<| < |r>|. In this way, the conver- 
gence of the Taylor expansion is guaranteed pro- 
vided that / is analytical everywhere except at 
the origin. Thus, the non-analyticity at the ori- 
gin explains why B functions as well as all the 
other commonly occurring exponentially decaying 
functions have addition theorems with a two-range 
form of the type of the Laplace expansion (|l.2| ) of 
the Coulomb potential. 

The crucial step, which ultimately makes the 
Taylor expansion method practically useful, is 
the expansion of the translation operator e'"<'^> 
in terms of differential operators which are irre- 
ducible spherical tensors. The starting point is an 
expansion of e'^ '' with a, b G in terms of mod- 
ified Bessel functions and Legendre polynomials 



(p. 108 of |& 



^ab 



„ab cos 9 



/ vr \i/2 



^ {21 + 1) /,+i/2(a6) P,(cos 9) . (3.3) 



Next, the series expansion for the modified Bessel 



function 1^^1/2 (p. 66 of 1 94 1) is inserted into ( |3.3| ), 



and the Legendre polynomials are replaced by 
spherical harmonics. This yields the following ex- 
pansion of e**'^ in terms of regular solid harmonics, 
which are irreducible spherical tensors of rank i, 
and even powers of the vectors a and b, which are 
irreducible spherical tensors of rank zero: 



i=0 m= 



A:=0 



2£+2fe^!(l/2; 



(3.4) 



This relationship is obtained from ( |3.3| ) by rear- 
ranging the Cartesian components of the vectors 
a and b. Accordingly, it holds for essentially arbi- 
trary vectors a and b and we can choose a = r< 
and b = V>, yielding 



E [^r(r<)]* 3^r(v= 



£=0 m=-t 
oo 



„2k y2fc 



< 



> 



to 2^+''fe!(l/2)m+i 



(3.5) 



It seems that (|3.5| ) was first published by San- 
tos (Eq. (A. 6) of [|105|| ). Santos also emphasized 
that this expansion should be useful for the deriva- 
tion of addition theorems, but apparently did not 
use it for that purpose. 

The only non-standard quantity in (^^) is the 
differential operator 3^^(V). Therefore, we will 
now discuss those properties of 3^7^(V) that are 
needed for the derivation of addition theorems for 
B functions. Other properties as well as numer- 
ous applications - mainly in the context of mul- 
ticenter integrals - are described in the literature 
|,|1,|9|,|1,|53,[ 



105Hl35|. 



The spherical tensor gradient operator y^{V) 
is an irreducible spherical tensor of rank (. (p. 312 
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of ||5J). Consequently, its application to a func- 
tion (j){r), which only depends on the distance r 
and therefore is an irreducible spherical tensor of 
rank zero, yields an irreducible spherical tensor of 
rank £ according to 



r dr / 



re 



(3.6) 



As discussed in Section IV of [112|, this relation- 
ship can be derived with the help of a theorem 
on differentiation which was published by Hobson 
already in the late 19th century |106| . Further de- 
tails can be found on pp. 124 - 129 of Hobson's 
book which was first published in 1931. 

If the spherical tensor gradient operator is ap- 
plied to a spherical tensor of nonzero rank, i.e., to 
a function that can be written as 



The abbreviations Al, A£i, Ai2, and a{i) are de- 
fined in (imi) - (|2l^ ). 

These two as well as analogous other expres- 
sions for 7£j^2 (^) '^^^ used for the derivation of 
addition theorems of B functions. However, in the 
case of B functions more convenient expressions 
are available. Consequently, the general expres- 
sions ( p.9p and ( |3.10| ) will actually not be used at 
all in this article. 

As mentioned before, 3^™(V) is obtained from 
the regular solid harmonic (r) by replacing the 
Cartesian components x, y, and z of r by the 
Cartesian components d/dx, d/dy, and d/dz of 
V. Consequently, 3^™(V) and 3^™(r) must obey 
the same coupling law. Hence, (^^) implies (Eq. 



(6.24) of PI ): 



F-^(r) = /,,(r)y-^(r/r). 



(3.7) 



the structure of the resulting expression can be un- 
derstood in terms of angular momentum coupling 
(Eq. (3.9) of pi): 



= ^ ^"^^ {imi + m2\£imi\£2m2) 
f=f ■ 



3^,7(V)F,7(r 



(2) (£m^ + m2\£imi\£2 



m2 



It follows from the summation limits ( |2.9D that the 
power £1 -\- £2 — £ = 2A£ of V is either zero or an 
even positive integer. 

Let us now assume that a spherical tensor 
F™^(r) and a radially symmetric function $^2(r) 



X lke2 (^) l^r'^"' (r/^) • (3-8) e^ist which satisfy 



For the radial functions J^-^^i^ in (^^), various rep- 
resentations could be derived, for example (Eqs. 
(3.29) and (4.24) of p3| ) 



3^™^(V)cI>,,(r). 



(3.12) 



Ae 



^ {-A£U-a{£)-l/2), 



q=0 



r dr 



El -I? 



Mr) 



(3. 



A^2 

E 

s=0 



{-A£2)s{A£i + l/2)s ,sa,~i,-2s~i 



2' r' 



If we apply 3^"^(V) to F;^'(r), then the two dif- 
ferential operators can be coupled according to 
(|1|). With the help of (jsj), we then obtain 
(Eq. (3.9) of m): 

3^,7(V)F-^(r) 

= ^^'^ {£mi+m2\£imi\£2m2)V^'^^^-^ 



1 d 
r dr 



J2+1 



(3.10) 



1 ^p.Ar] 

r dr 



mi+m2 / 



(3.13) 
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This relationsiiip is particularly well suited for 
B functions. In this case, ( 3.13| ) is more conve- 
nient than other, more general expressions for the 
product {'V)F^^ (r) which can for instance be 
found in articles by Santos |105], Dayman |107|, 



Stuart I 111 I, Niukkanen |115], Weniger and Stein- 
born pi, and Rashid [I2C 



4 B Functions 

In this Section, those mathematical properties of 
B functions, that are relevant for the derivation of 
an addition theorem via three-dimensional Taylor 
expansion, will be reviewed. More complete treat- 
ments can be found in p|, [lT^ , |l2^ , p^ |l40|. 

If Ky{z) is a modified Bessel function of the 
second kind (p. 66 of the reduced Bessel 

function is defined by (Eqs. (3.1) and (3.2) of 

M) 



(4.1) 



If the order of a reduced Bessel function is 
half-integral, v = n+1/2 with n € No, the reduced 
Bessel function can be written as an exponential 
multiplied a terminating confluent hyper geometric 
series iFi (Eq. (3.7) of [|l3^ 



K+l/2{z) 

= 2"(l/2)„e-^Fi(-n;-2n;2z) 



(4.2) 



The polynomial part in ([4. 2]) was also treated in- 



dependently in the mathematical literature | 142 |, 
where the notation 



Qn{z) 



■ 



(4.3) 



is used. Together with some other, closely related 
polynomials, the G„(z) are called Bessel polyno- 
mials. They are applied in such diverse fields as 
number theory, statistics, and the analysis of com- 



plex electrical networks |142|. 

The so-called B function was introduced by 
Filter and Steinborn as an anisotropic general- 
ization of the reduced Bessel function (Eq. (2.14) 
of |9|), 

= [2"+^(n + 0!]-^A;„-i/2(ar)3^r(«r), (4.4) 



where a > and n G Z. Because of the facto- 
rial {n + l)\ in the denominator, B functions are 
defined in the sense of classical analysis only if 
n + I > Q holds. However, the definition of a S 
function remains meaningful even for n + d < 0. 
If r 7^ 0, -B™„_££(r) with n G N is zero, but for 
r = 0, its value is oo/oo and therefore undefined. 
In fact, such a B function can be interpreted as 
a derivative of the three-dimensional Dirac delta 
function (Eq. (6.20) of [pTp : 



(2^- l)!!47r 



a 



£+3 



X [l-a-^V^] dj'ir), nGN. (4.5) 
The spherical delta function 5"^ is defined by 



i-lY 



(2£- 1)!! 



3;r(V)<5(r). (4.6) 



If follows from ( pT| ) - (|J) that B functions 
are relatively complicated mathematical objects. 
Nevertheless, all the commonly occurring expo- 
nentially decaying functions can be expressed as 
simple finite sums of B functions. For example, if 

x:^,(a,r) = (ar)"-ie— yr(^,0) (4-7) 

is an (unnormalized) Slater-type function with 
a > and n > ^ -Fl, then (Eqs. (3.3) and (3.4) 
of |9l) 



n-e 



l-p 



P=Pn 



{2p - n + i)\{2n - 2i - 2p)\] 



BZ{a,v), (4. 



where Pmin — 

{{n - i)/2\. Here, {x} stands for 
the integral part of x, i.e., for the largest integer 
V satisfying v < x. 

As is well known, the bound-state eigenfunc- 
tions of a hydrogenlike ion with nuclear charge Z 
can be expressed as follows. 



W^n":.(^,r) 



X e-^^/" LfA,l\{2Zr/n) yr{2Zv/n) , (4.9) 





"(n-£- 1)!' 


3/2 




2n{n + t)\ 
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where L^j^^^^\ is a generahzed Laguerre polyno- 



mial. Then (Eq. (3.16) of |3|), 

n J 

n{n + (. 



2^+1 



(2^ + 1)!! 

n-e-i 



1/2 



2{n-£-l)\ 
i-n + i + l)t{n + e + l)t 



t=o 



tl{i + 3/2)t 



(4.10) 



The bound-state hydrogenic eigenfunctions are or- 
thonormal, but not complete in the Hilbert space 
L^(M^). This incompleteness is sometimes over- 
looked (see for example the discussion in |143| ). 

The following set of functions, which was in- 
troduced by Hylleraas |144] and by Shull and 
Lowdin | 145| , 146 1, is complete and orthonormal 
in the Hilbert space L'^' 



A™,(a,R) = (2a) 



3/2 



IV. 



(n + e + i) 



e-"^4'5-!(2«03^r(2ar). (4.11) 



Then (Eq. (3.18) of p7| ), 

A™,(a,r) = (2a)3/2 



X 2 



(2n + 1) 
(2^ + 3)!! 

' E - 

t=0 



(n + ^ + 1)! 



1/2 



n-l-\)\ 
n + ^+l)t(n + £ + 2)t 



t!(^ + 5/2)t 



xS- ,(a,r) 



(4.12) 



Closely related to the bound-state hydrogenic 
eigenfunctions is the following set of functions 



which were already used in 1928 by Hylleraas |148| 
and which are commonly called Coulomb Sturmi- 
ans or simply Sturmians |]l4 



*-,(a,r) = (2a) 



3/2 



n 



2n(n + i)l 



3/2 



X e-"'- Ll'X[i2ar) 3^r(2ar) . (4.13) 



Comparison of ( |4.9D and ( |4.13D yields 

f-,(Z/n,r) = W:re{Z,r). (4.14) 



The fact, that bound-state hydrogenic eigenfunc- 
tions and Sturmians have the same normalization 
constant, is actually by no means obvious. The 
W^^{Z, r) are normalized with respect to norm of 
the Hilbert space L^{R^), whereas the ^'™^(a,r) 
are normalized with respect to the norm of the 
Sobolev space W^iR^) m. Further details can 
be found in []53| . 

In ( 4.1C| ), we only have to replace Z/n by a to 
obtain the following representation of Sturmians 
in terms of B functions (Eq. (4.19) of p3|): 



M/™,(a,r) = (2a)='/^ 



(2^ + 1)!! 

n-l-l 

t=o 



n{n + iy. 



1/2 



2(n-£ - 1)! 
{-n + i + l)t{n + i + l) 



tKi + 3/2)t 



xBr+,/a,r). 



(4.15) 



The hnear combinations (U), ( [4.10| ), ( [1.12D , and 
( [4.15| ) imply that mathematical results for B func- 
tions can be translated immediately to analogous 
results for Slater-type functions, bound-state hy- 
drogenic eigenfunctions. Lambda functions, and 
Sturmians. 

In view of the fact that B functions have a rel- 
atively complicated mathematical structure, it is 
by no means trivial that the linear combinations 
(PI) , (l4lc| ), (|1|), and (|1|) could be derived at 



all. However, these relationships as well as other 
advantageous properties of B functions can be ex- 
plained via their Fourier transform, which is is 
of exceptional simplicity among exponentially de- 



clining functions (Eq. (7.1-6) of [137] or Eq. (3.7) 

of inr 



^™,(p) = (2^) 



-3/2 / e-*P'-S™^(r)d^r 



(2/vr: 



,1/2 



a 



2n+e-l 



[a2 + p2]n+£+l 



yn-ip)- (4.16) 



The h near combinations (U), ( [4.10| ), (|]T2|), and 
( 4.15| ) imply that the Fourier transforms of Slater- 
type functions, of the bound-state hydrogenic 
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eigenf unctions, of Lambda functions, and of Stur- 
mians can be expressed as simple linear combina- 
tions of Fourier transforms of B functions. 

If we want to derive an addition theorem for 
some function via the translation operator in its 
spherical form (|3.5| ) , then the application of higher 
powers of the Laplacian and of the spherical ten- 
sor gradient operator to this function must lead to 
expressions of manageable complexity. On the ba- 
sis of the Fourier transform ( 4.16| ), it can be seen 
that these expressions are indeed very simple in 
the case of B functions. 

The differential operator 1 — V^/a^, which is 
typical of the modified Helmholtz equation, acts 
on B functions as a ladder operator for the order 
n (Eq. (5.6) of p3[): 

[l-Vya']BJ^Ar) 



(r) . (4.17) 



This follows at once from the Fourier transform 
(|4.16| ) in combination with the fact that in mo- 
mentum space the differential operator [1 — V^/a^] 
is transformed into the multiplicative operator 

With the help of this relationship, the applica- 
tion of the higher powers of the Laplacian can 
be expressed quite easily in closed form. Combi- 



nation of (4.17) with the binomial theorem yields 
(Eq. (5.7) of 11131): 



^—2u \-i2v Tjm /„\ 



i=0 



(4.18) 



The remarkably simple Fourier transform 
( 4.16| ) also explains why the application of the 
spherical tensor gradient operator to a -B function 
produces particularly compact expressions. For 
example, if we combine 

yriV)e^'P-'' = y^{±ip)e^'P-'' (4.19) 

with the Fourier transform ( |4.16| ), we find that a 
nonscalar B function can be obtained by applying 
the corresponding spherical tensor gradient oper- 
ator to a scalar B function (Eq. (4.12) of [lll2| ): 



This relationship can also be obtained in a 
straightforward way with the help of ( |3.6D . The 



differential operator (r~^d/dr) in (3^) is a sim- 
ple operator for modified Bessel functions, but not 
for the other exponentially decaying functions con- 
sidered in this article. It acts as a shift operator 
for reduced Bessel functions according to 



1 H \ 

k{z) = (-irk.miz). (4.21) 



Similarly, the application of the spherical ten- 
sor gradient operator to a nonscalar B function 
with ^ > yields a simple linear combination of 
B functions (Eq. (6.25) of |ll3| ): 

^ /AA 



This relationship follows at once from ( 3.13| ), 
(|1|), and ( pO|) . 



5 An Addition Theorem for B 
Functions 

We now want to derive an addition theorem for B 
functions via three-dimensional Taylor expansion. 
For that purpose, we apply the translation opera- 
tor in its spherical form ( |3.5D to a i? function: 



^"£(«,r< +r>) 

OO l\ 

^ rl^V 

q=0 



3^r(v>) 



S 2^(1/2)..... ^"":-^"'->^-^^-^^ 
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Now, we want to apply the powers of the Lapla- kind according to (p. 66 of 
cian V> to the B functions in (^.ip. This can be 



done with the help of ( [4.18| ), yielding 
S™,(a,r<+r>) 



3Fi(z. + l;zV4) 



r(z^ + i) 

{z/2r 



(5.6) 



E [yr: 

oo 



(ar<)29 



q=0 



2^1+2,^,(1/2)^^ 



By combining (p.5|) and (^■6p we obtain 



E(-l)*(!)^"-MKr>). (5.2) X E(-l)'0^-.Kr>) 



It follows from the definition of a B function ac- 
cording to (gj) that 5™^(r>) with r> 7^ is 
zero for n < —i. Thus, the summation limit of 
the innermost summation in ( ^.2| ) is not q but 
min(g,n + i). Hence, if we introduce a new sum- 
mation variable s = q — t, we obtain for the two 
innermost summations in (B.2|): 



(ar<) 



2q 



q=0 



(1/2)., 



+q+l 



X 



^(-1)* Qi?-*/"'^; 

(ar<)2 



^ 2^i+2.g!(l/2),,+,+i 

min(g,n+^) 

X 



5^ (-l)*(^Jji?-,,,(a,r>) (5.3) 

"+^ 



s=0 



(-l)^(ar</2) 



^ ^ t!(l 



2t 



^ t!(l/2),,+i+i 
oFi(^i+t + 3/2;(ar<)V4). (5.5) 



-1/2?- (-l)*2i/2-* , , , 

- E t!(ar<)^i-m/2 ^^i+m/2(«r<) 

xS- ,(«,r>). (5.7) 



Inserting ( |5.7D into (|5.2D yields: 

^^£(a,r< +r>) 

00 l\ 

= (2-)^/^E E i^r* 

^1=0 m=-£i 



E^(«^<)*"'^"'^'^^i+m/2(ar<) 



(5.8) 



Now, all that remains to be done is the applica- 
tion of the spherical tensor gradient operator to 
the B functions. This can be done with the help 
of (|422|) , yielding 



3^-HV>)i?r-M(r>) 

= (-a)^i ^ ^'^^ (i>2m^mx\hmx\i' 



M2 



E(-ir i^l') B-:r-i-s-tM^->) (5-9) 



The generalized hypergeometric series 0-^1 can be 
replaced by a modified Bessel function of the first 



The summation limit in (|5.9| ) is defined in 
( |2.13| ). By inserting (|5^ ) into (^), we finally 
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obtain an addition theorem for B functions: 
^™£(a,r< + r>) 

£1=0 m=-ei 



t=0 



X E(-i)M , pr+n-.-M.(r>)-(5-io) 

This expression is not yet completely satisfactory 
from a formal point of view. We have to take into 
account that the B functions in the innermost sum 
are zero ifn + £ — £2 — s — t < —i2- For that 
purpose, we introduce a new summation variable 
q = n + i — t. This yields: 



(27^)3/2 

(_2)n+£ 



E E (-1)'^ [Knr<) 

£1=0 m=-£i 



E 

g=0 



(-2)^ 



(n + . 



vn+£-£i-9-l/2 



X ^n+£+£i-q+l/2 



min(5,A£2) 

E (-1 



m) 



s=0 



A£2 



(5.11) 



This addition theorem is identical to the ad- 
dition theorem (4.27) of which was derived 
with the help of techniques that closely resemble 
the zeta function method of Barnett and Coul- 



son |31|. Other addition theorems for B functions 
and for reduced Bessel functions can be found in 
Sections IV and V of 



6 Summary and Conclusions 

Addition theorems are expansions of a given func- 
tion /(ribr') with r, r' G in terms of other func- 
tions that only depend on either r or r'. As dis- 
cussed in Section ||, many different techniques for 
the construction of addition theorems are known. 
Examples are the zeta function method of Bar- 
nett and Coulson |31|, Lowdin's alpha function 
method [p3] , or the Fourier transform method sug- 



gested independently by Ruedenberg |51] and Sil- 
verstone [52|. 

Addition theorems can also be defined as 
three-dimensional Taylor expansions according to 
( 11.31 ). In this approach, the translation opera- 
tor e"" generates f{v + r') by doing a three- 
dimensional Taylor expansion of / around r. In 
the context of atomic and molecular calculations, 
one is normally interested in addition theorems of 
irreducible spherical tensors of the type of (|1 . 1| ) . 
In such a case, the straightforward application 
of the translation operator in its Cartesian form 
gX d/dx^y d/dy^z d/dz -^q^^^j iga,d to enormous tech- 
nical problems. Consequently, it was generally 
believed that three-dimensional Taylor expansions 
constitute essentially a formal solution, but can- 
not be used for the actual derivation of addition 
theorems of irreducible spherical tensors. How- 
ever, the Taylor expansion method is indeed a 
practically useful tool if the translation operator 
is expanded in terms of differential operators that 
are irreducible spherical tensors according to ( |3.5| ) . 
This was demonstrated in ||5^ by constructing the 
Laplace expansion of the Coulomb potential as 
well as the addition theorems of the regular and 
irregular solid harmonics and of the Yukawa po- 
tential. 

Since the addition theorems mentioned above 
are all comparatively simple, their successful con- 
struction does not necessarily guarantee that the 
Taylor expansion method is also practically use- 
ful if more complicated addition theorems are to 
be constructed. Consequently, the topic of the 
present article is the construction of addition the- 
orems for exponentially decaying functions which 
are notoriously complicated. 
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Slater-type functions have the simplest struc- 
ture of all the commonly occurring exponentially 
decaying functions. Consequently it looks like an 
obvious idea to concentrate on addition theorems 
of Slater- type functions, and to express the ad- 
dition theorems for other exponentially decaying 
functions in terms of addition theorems for Slater- 
type functions. 

However, the addition theorems for Slater- 
type functions are notoriously difficult. More- 
over, Slater-type functions have a simple analyt- 
ical structure only in the coordinate representa- 
tion. In momentum space, which is probably more 
important in the case of multicenter problems, 
the simplest exponentially decaying functions are 
the so-called B functions. Their extremely simple 
Fourier transform ( [4.16| ) as well as some other ad- 
vantageous mathematical properties give B func- 
tions a unique position among exponentially de- 
caying functions. Moreover, as discussed in Sec- 
tion all the commonly occurring exponentially 
decaying functions can be expressed as simple fi- 
nite sums of B functions. This implies that the ad- 
dition theorems of Slater-type functions, of bound- 
state hydrogenic eigenfunctions, and of other func- 
tions based on Laguerre polynomials can be ex- 
pressed in terms of the addition theorems of B 
functions. 

So, this article focuses on the derivation of an 
addition theorem for B functions by doing a three- 
dimensional Taylor expansion in its spherical form 
(3.5). The relative ease, with which the addition 
theorem ( 5.1l| ) can be constructed, and its com- 
paratively compact structure indicate once more 
that B functions have in the context of multicen- 
ter problems considerable advantages over other 
exponentially decaying functions. 
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